Switchable lasing in coupled multimode microcavities 
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We propose the new concept of a switchable multimode microlaser. As a generic, realistic model of a multi- 
mode microresonator a system of two coupled defects in a two-dimensional photonic crystal is considered. We 
demonstrate theoretically that lasing of the cavity into one selected resonator mode can be caused by injecting an 
appropriate optical pulse at the onset of laser action (injection seeding). Temporal mode-to-mode switching by 
re- seeding the cavity after a short cool-down period is demonstrated by direct numerical solution. A qualitative 
analytical explanation of the mode switching in terms of the laser bistability is presented. 
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Microlasers with cavity sizes comparable to the radiation 
wavelength A are very promising from both fundamental and 
application points of view |1] for use as integrated coher- 
ent light sources. Making microlasers capable of multiple- 
wavelength emission contributes even more towards miniatur- 
ization of optical components, and also provides an additional 
degree of freedom in light control. The common approach to- 
wards microlaser tunability is in essence modification of the 
optical properties of a single-mode cavity by thermal mi- 
cromechanical [3] or electrooptical means [4]. 

In this Letter we propose, in contrast, the new concept of 
a switchable microlaser, comprised of a multimode laser mi- 
croresonator, where lasing can be switched on demand to any 
of its eigenmodes. While the broad gain profile of semicon- 
ductor and dye lasers provides little discrimination between 
neighboring modes and, thus, leads to mode hopping and mul- 
tistable mode dynamics, we show that a definite resonator 
mode can be selected for lasing by injection seeding [5, d], 
i.e. by injecting an appropriate pulse before and during the 
onset of lasing, such that the stimulated emission builds up 
in a chosen mode from this seeding field rather than from 
the random noise present in the system due to quantum fluc- 
tuations and spontaneous emission. We investigate the time 
needed for switching between different lasing modes and an- 
alyze how this process is influenced by noise in competition 
with the seeding signal. 

We generalize the semiclassical multimode laser model of 
Ref. ITI] to study the dynamics of individual modes in a two- 
mode laser, where an external field is present in the cavity 
due to the injection seeding. Semiclassical theories of optical 
amplifiers [5] have been successfully employed to study the 
mode dynamics in multimode lasers [7, 8], especially in the 
context of laser instabilities [9] . An example of a two-mode 
microcavity, which we have chosen to study numerically, is 
depicted in Fig. [T]a). The microcavity is based on two cou- 
pled defects in a two-dimensional (2D) photonic crystal (PhC) 
made of dielectric nanopillars [10]. Pillar PhCs are practically 
feasible by state-of-the-art fabrication technology 1 1 1]. 

The system of two identical coupled single-mode cavities 
supports two modes, namely the bonding (symmetric) and the 
antibonding (antisymmetric) mode, characterized by spatial 



field distributions ^^1,2(1*) and frequencies uji^2 = ^0 T 
respectively. Here Aco is the mode detuning from the fre- 
quency of the single-cavity resonance, coq. For weak mode 
overlap the spatial intensity profiles of two modes nearly 
coincide, |^^i(r)| ^ |^2(^)| • We assume that the cavi- 
ties contain a laser medium with the gain centered at fre- 
quency uja = ujQ -\- homogeneously broadened to have a 
line width of AcOa > Aco, where S is the detuning of the gain 
profile from the cavity frequency coq. See Fig.[T]b) for the def- 
inition of the frequency parameters. For the numerical calcu- 
lations below we take a quadratic array of 12 x 12 nanopillars 
with dielectric constant e = 9 in air (e = 1). The lattice pe- 
riod is a = 500 nm and the pillar radius r = 0.2a. Then both 
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Figure 1: (Color online) a) Schematic illustration of the resonant 
system under study. Two coupled defects in a periodic 2D PhC lattice 
contain a gain medium. Two emitters shown are used to produce 
the seeding signal delivered to the defect sites by two waveguides. 
The spatial field distribution of the seeding signal is schematically 
depicted as a gray scale on the left, b) Frequency parameters of the 
two-mode laser, see text. 
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defect modes have frequencies near a/X^ 0.395 [Fig.[T]b)]. 

To obtain an initial understanding of the seeding-induced 
mode switching, the semiclassical Max well-B loch equations 
[5], may be simplyfied using the rotating-wave and the slowly 
varying envelope approximations. All the spatial dependen- 
cies of the electric field and atomic polarization can then be 
represented in the basis of the two cavity modes, such that 
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atomic polarization can be eliminated adiabatically ItIS]. For 
class-A lasers, where the radiative decay rate 7^ = Aa;a/2, 
the non-radiative decay rate 7||, and the cavity mode decay 
rates hij are related as 7^ ^ 7|| ^ f^j^ the following system 
of equations for the slowly varying envelopes Ej {t) in the two 
modes is obtained i?!] j = 1, 2, i j), 
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In Eqs. ([B the terms linear in Ej (t) describe stimulated emis- 
sion driving and are controlled by the light-matter coupling 
g y/27ru;od'^ / h, with d the dipole moment of the atomic 
transition, by the pumping rates projected onto the modes 
j = 1,2, Rj = J^Uj{T)uj{T)R{T)dT, and by the cav- 
ity mode decay rates tZj. The coefficients Cj = Im/3~^, 
with (3i^2 = S ± Aco — zAcJa/2, account for the different 
mode-to-gain couplings due to asymmetrical detuning of the 
atomic transition with respect to the resonator lines. The 
terms cubic in Ej (t) describe field saturation above the las- 
ing threshold, where r] = d^ /2^\\h^ and the overlap integrals 
a^i = u* {r)uj {r)ul {r)ui (r)(ir are taken over the regions 
G containing the gain medium. Since \ui{r)\'^ ^ |^2(^)|^ we 
can assume aj'j = aj- = a, R 

The cross-saturation terms, with M 



R2 = R and ni = ^2 = k,. 
y =/3-i + and 
Xi,2 = — nil / — ijw), depend in an assymetrical way 

on the mode indeces j. However, this asymmetry remains 
small unless Auj <C Acj^. 

The inhomogeneous terms Fj (t) originate from the exter- 
nal injection seeding field and from a noise field account- 
ing for spontaneous emission 0]. For vanishing functions 
Fj, Eqs. ([B would take the form of the standard two-mode 
competition equations |5] with mode coupling constant C 
slightly exceeding unity. This corresponds to bistable lasing 
lll2n and to mode hopping in the presence of stochastic noise 
in the system [13]. If both an external seeding field ^*(r, t) 
and a stochastic noise field £^^(r, t) are present in the cavity, 
f (r, t) = S'{r, t) + f ^(r, t), one obtains. 
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= ™)+Ff(t). 



The time integration in Eq. Q is the averaging over a time 
interval larger than 1/Acj. The function F{t) is determined 
by the temporal dependence of f*(r, t). The coefficients 
Fj and Fj^{t) are determined by the spatial overlap of each 
mode with the seeding and noise fields, respectively. We con- 
sider the situation when the seeding prevails over the noise, 
i.e., FjF{t) Ef{^)^ before and during the onset of lasing. 



After the onset the Ej become so large that the terms Fj have 
no effect anymore. In this situation the evolution of the res- 
onator will be determined by the ratio of and F| . 

We analyze the influence of the balance between Fl and 
on the final lasing state by numerical solution of Eq. ([T]) and 
plotting the phase trajectories of the temporal resonator state 
evolution in the (|Fi|^ , |F2|^) plane (Fig. [21). The numeri- 
cal values of the coefficients in Eq. ([T]) are calculated for the 
model PhC structure of Fig. [T] a). As seen in Fig.[2l the lasing 
state first reaches overall intensity saturation ( |Fi | ^ + | F2 1 ^ = 
El) and then drifts towards one of the stable fixed points 
corresponding to single-mode lasing (either \Ei\^ = E^ or 
1^2!^ = Eg). The drift occurs after the sharp bend seen in 
each of the phase trajectories. If the mode coupling constant 
is C = 1, the drift becomes infinitely slow, and the (1,0)-(0,1) 
line in Fig. [2] turns into a line of fixed points. For our case, 
where C only slightly exceeds unity, the drift happens on a 
longer time scale than the initial overall intensity growth, and 
the intermode beats decay fast after the lasing onset. In the 
case of symmetric detuning of the cavity modes with respect 
to the gain frequency (S = 0), single-mode lasing is achieved 
into that mode whose spatial overlap with the seeding field. 
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Figure 2: (Color online) Cavity phase diagrams for a lasing system 
governed by Eq. flj for F|F(t) > F^{t) for (a) symmetric (6 = 0) 
and (b) non-symmetric (6 = — 0.25Aa;) values of the mode fre- 
quencies with respect to the central gain frequency to a ((^a < a;o, as 
shown in Fig.[T](b)). The dots denote the stable cavity states, and the 
curves represent the phase trajectories for their temporal evolution 
for different ratios Ff : F2 in the direction of the arrows. 
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Fj, is largest [Fig. [2] a)]. The asymmetry of the modes with 
respect to gain (6 7^ 0) shifts the turning point towards one of 
the modes, but if the seeding is chosen in a way that the spatial 
overlaps in Eq. (|2]) result in Ff ^ or Ff <C Fl , each of 
the modes can nonetheless be selected for lasing [Fig.[2lb)]. 

To refine the predictions of this simple theory of bistable 
lasing we have modeled the lasing action in coupled PhC de- 
fects with a realistic injection mechanism (Fig.^ using the 
finite-difference time-domain (FDTD) method 10, • The 
defect modes are located inside the PhC band gap. Both de- 
fects are filled with an active medium whose population dy- 
namics are described at each space point by the rate equa- 
tions of a four-level laser with an external pumping rate Wp. 
In defining the model and its parameter values we follow in 
detail Refs. In particular, the non-radiative transi- 

tion times of this model are taken so as to achieve popula- 
tion inversion, i.e. rs2 — no <^ ^21, with rai = rio = 
1 X 10~^^ s, r2i = 3 X 10~^^ s, and the total level population 
is A^totai = 10^^ per unit cell [16]. The Maxwell equations, 
supplemented by the usual equation of motion for the polar- 
ization density in the medium and by the laser rate equations 
il4 |15L Il6l I17II . are solved for the geometry of Fig. [T] a) in 
TM polarization, where E(r, t) = Ey{x^ t) y. The seed- 
ing signal is excited by two emitters (linear groups of dipoles) 
engineered on the same chip as the PhC and is transmitted to 
the defects through waveguides in the PhC (see Fig.[T]). Each 
of the emitters generates a single short Gaussian pulse with 
carrier frequency co at or near coa- The calculations have been 
performed using different, fixed values of the half- width du- 
ration at in the range between = 5 x 10~^^ and 10~^^ s. 
The relative phase of the fields in these pulses is chosen or tt. 
As expected, such seeding patterns almost exclusively excite 
the bonding and antibonding mode, respectively. Technically, 
the seeding dipoles are realized as pointlike oscillating current 
sources in the Maxwell equations Q 8] . Similarly, the sponta- 
neous emission (isl [l^ is modeled as an ensemble of 
point current sources, randomly placed in space, with tempo- 
rally (^-correlated Langevin noise [19]. For the FDTD compu- 
tations the computational domain 13a x 13a with perfectly 
matched layer (PML) boundary conditions was discretized 
with an a/ 16 mesh and a time step of dt = 6 x 10~^^ s re- 
quired by numerical stability, see 112 ill for details. 

Neglecting noise at first, we find that, if the laser amplifica- 
tion line spectrally covers both resonant modes and provides 
a comparable effective gain for each of them (i.e. |^| < Acj), 
a seeding signal of the type described above can individu- 
ally select any of the two modes when applied (with arbitrary 
strength) during the onset of lasing, i.e. during the exponen- 
tial growth after switching on the pumping. The steady- state 
lasing is then nearly single-mode, and the dominant mode is 
the one whose symmetry matches that of the seeding signal 
(Fig. [3]). This is in agreement with the semianalytical the- 
ory described above. The electric field maps in Fig. [3] (bot- 
tom) show that the spatial field distribitions of the two lasing 
modes remain nearly unaffected by the side-coupled seeding 
waveguides. In the present setup the laser output is primarily 
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Figure 3: (Color online) Illustration of mode-to-mode switching. 
Top-, time dependence of the normalized pumping power Wp, seed- 
ing signal 2y^ and cavity field Ey. Bottom', spatial distribution of 
the electric field amplitude Ey(x, z) (grey scale) in the steady-state 
lasing regime after initial seeding (1) and each re-seeding (2, 3). 
The time instants (l)-(3) are marked in the top panel by dashed 
lines. This calculation was done for a typical pumping rate Wp — 
1 x 10^^ s~^ and a seeding pulse duration at = 1.2 x 10~^^ s. 



delivered through the seeding waveguides. A detailed study 
of the microlaser radiation characteristics will be given else- 
where. Moreover, we have shown that successive switchings 
from one mode to another are possible. To achieve this, the 
pump is first turned off to allow the currently lasing mode 
to decay. After a certain cool-down time Tc, the pumping 
is turned back on, and the cavity is re-seeded for the other 
mode in the same way as the initial seeding occurred. Fig. [3] 
demonstrates such a switching sequence from the bonding to 
the antibonding mode and back. The pump was kept on for a 
time needed to achieve quasi-continuous lasing, as confirmed 
by the temporal dependence of the cavity field. Looking at 
the spatial field distribution in this regime [panels (l)-(3) in 
Fig. m, we make sure that the mode switching occurs in the 
desired order. 

The switching time between two lasing modes is primar- 
ily controlled by the minimum cool-down time Tc^min^ i-^- by 
the relatively slow mode decay time, r = 1/ and hence by 
the cavity Q-factor Q = coj/k,, not by the fast lasing onset 
after re-seeding, Ts ~ Tc^min- Since the mode decays expo- 
nentially with time from its steady- state lasing amplitude Aq 
and the re- seeding signal must be strong enough to override 
the residue of the decaying initial mode, Tc^min decreases log- 
arithmically with increasing re- seeding pulse power (i.e. with 
its time-integrated intensity) or, for a Gaussian pulse, with the 
seeding amplitude S, Tc^min = ^) ^'^{S/Aq). Here, Cj 

is a dimensionless factor describing the coupling of the cav- 
ity mode j = 1, 2 to the laser transition. The corresponding 
numerical FDTD results are shown in Fig. |4l For realistic 
parameter values (Fig.|4l) switching times of a few tens of pi- 
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Figure 4: (Color online) The minimum switching time rc,min versus 
the re- seeding amplitude, normalized to the saturated lasing mode 
amplitude. The seeding pulse duration is at = 0.12 ps. The cavity 
mode decay time is r = 0.7 ps. The longer rc,min for switching 
from the bonding to the antibonding mode is due to the non-zero 
gain detuning S towards the bonding mode frequency [see Fig.[T]b)]. 



coseconds can be realized for re-seeding amplitudes 5* as low 
as 0.001 per cent of the lasing mode amplitude Aq, two to 
three orders of magnitude faster than in previous devices [22] . 

When noise is present in the cavity, the FDTD calculations 
show that controlled switching is preserved as long as the in- 
tegrated seeding power exceedes the noise power integrated 
over the lasing onset time. Otherwise, noise begins to dom- 
inate the lasing spectrum formation. In this case, the mode 
whose spatial overlap with the noise field at the onset of las- 
ing is larger wins the competition. The same happens when 
the seeding pulse does not match the onset of lasing in time 
and its field residue has smaller amplitude in comparison to 
the noise field. These numerical observations are fully con- 
sistent with the analysis of the competition equations ([B with 
the functions Fj{t) given by Q. The switching is effective 
only if the influence of the seeding prevails in the resonator at 
the period of time when laser radiation starts to build up. 

To summarize, the concept of switchable (rather than tun- 
able) lasing in microstructures has been introduced. Instead 
of externally changing the parameters of a single-mode cavity, 
an inherently multimode cavity is used, and one of the modes 
is deliberately made to be dominant for lasing by means of 
injection seeding. This offers the possibility of all-optical fre- 
quency selection and switching in microlasers with particu- 
larly low switching times. As an example, we have inves- 
tigated the mode switching in a system of two coupled de- 
fects in a 2D PhC lattice. For realistically chosen parameters 
a mode-to-mode switching on the picosecond scale has been 
numerically demonstrated. The results are consistent with a 
qualitative semianalytical model. It shows that a resonator 
supporting modes with similar spatial intensity profile tends 
towards bistability, which is the underlying physical mecha- 
nism of switchable lasing. The proposed concept is not lim- 
ited to the model considered, but is expected to work in any 



resonator featuring bi- or multistability. Any coupled cavity 
based system would be a good candidate for the effects pre- 
dicted. For instance, we have observed four-mode switching 
in coupled nanopillar waveguides of both periodic and non- 
periodic longitudinal geometry [181 . 
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